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ABSTRACT 



Context. We study the planet populations in the discovery window of the CoRoT-space-telescope scheduled for launch on December 
27th. We base the prediction on 'first principles' calculations of planet formation in the framework of the planetesimal hypothesis. 
Aims. To provide a-priori planetary initial mass functions for confrontation with the CoRoT-planet discoveries in the entire range of 
sensitivity of the CoRoT instrument, i.e. for all giant planets and down to terrestrial planet masses. 

Methods. We construct a comprehensive set of static complete-equilibrium core-envelope protoplanets with detailed equations of 
state and opacity and radiative transfer by convection and radiation. Protoplanets are calculated for host-star masses of 0.8 to 2 solar 
masses and orbital periods of 1 to 16 days. We subsequently check the stability of the planetary population by a series of methods. 
Results. We find the static planetary populations to be stable and thus a plausible ensemble to predict the planetary IMF for orbital 
periods in the specified range. 

Conclusions. We predict bimodal planetary initial mass functions with shapes depending on orbital period. The two main maxima 
are around a Jupiter mass and about 50 earth masses. We predict an abundant population of Hot Neptunes and a large population of 
planets that fill the solar-system gap of planetary masses between Neptune and Saturn. 

Key words. planets:formation - exoplanets - planets: mass - solar system:formation - initial mass function 



1. Introduction 

Planets form in circumstellar disks which are a necessary step 
in star-formation. Presently two mechanisms are in discus- 
sion: formation via a d isk-instability and via planetesimals, see 
IWuchterl et al.l d2000t) . Uncertainties about the actual nebula 
conditions at the era of planet formation and an incomplete 
understanding of relevant physical processes (e.g. of planetes- 
imal formation and of the various types of orbital migration) 
limit the predictive power of conventional formation theories. 
In preparation for t he CoRoT-transit-search for exoplanets, e.g. 
Bagl in et ail d2002l) we have develo p ed a new synoptic approac h 
to planet form ation, dPecnikl d2003l) : iPecnik & Wuchterl((l2005l) : 
Broeg d2006j)). 



2. Synopsis of protoplanets 

Broeg] d2006l) calculated all planets and protoplanets in hydro- 
static and thermal equilibrium for stellar masses of 0.4, 0.8 , 
1 and 2M G Q and orbital periods of 1, 4, 16 and 64 d (0.02, 
0.049, 0.124 and 0.313 AU for the solar case). The ranges are 
chosen to cover CoRoT's discovery space. Planetary models 
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1 The stellar types that can be observed with high photometric pre- 
cision and are likely to be found in CoRoT's fields. The period range 
results from the 150d continuous observation periods and the require- 
ment to observe three transits to make a safe detection. 



are spherically symmetric, consist of a rigid, constant density 
core and a gaseous envelope in hydrostatic and thermal equi- 
librium. Planetesimal accretion provides an energy input at the 
core surface that is transferred through the envelope by radia- 
tion and convection. Detailed equations of state and opacities 
for roughly solar composition are used with 'interstellar' dust 
opacities. Orbital distance enters via the thermal and tidal ef- 
fects of the star. Assuming that no core size is favoured and that 
gravitationally stable nebulae provide the pressures for embed- 
ding the planets, statistical properties of the planet population 
are provided, in particular planetary mass spectra. Here we use 
the Co RoT Mark 1 vl.l ensemble of planetary models, Broeg 
(2006), Mark 1 for s hort. Mass spectra and solution manifolds 
are av ailable online, iBroeg. C.1 (|2006), see Broeg & Wuchterll 
(2006) for an application to the large-core Pegasi-planet HD 149 
026 b. Here we apply the method to theoretically predict what 
planets are in CoRoT's search space. We restrict the discussion 
to close-in planets with orbital periods of < 16 d for two reasons: 
(1) they are in CoRoT's discovery window for the first run with 
a foreseen continuous observing period of 30-60 d and hence 
with a completeness limit for planets detected by at least three 
transits at about 15 days, (2) there is a qualita t ive ch ange in the 
planet-population for periods of ^ 16 d. lBroegl ((2006) and planet 
formation being dominated by hyd rostatic processes for smaller 
periods, IBroeg & Wuchterll d2006l) . 
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Fig. 1. Theoretical planetary initial mass functions calculated from planet formation theory for a typical CoRoT-field. Results are 
shown for planetary orbital periods of 1, 4, and 16 days, from left to right. The relative frequency is plotted as function of lg mass 
in kg. Vertical red lines mark the Earth and Jupiter masses. ~ 10 6 planetary models in total. Structures of width < 0.3 dex have to 
be taken with care, because of undersampling in spectral type due to the unexpected richness of the mass-spectra. ('M2+M4+M6' 
designates planetary core-accretion and is not related to the stellar population). 



3. Planetary Mass Spectra for a typical CoRoT field 

Our calculations of the planet population show a strong depen- 
dence of the planetary properties on stellar mass, lBroeg| (l2006). 
We therefore account for the types of stars CoRoT will observe. 
Given the Mark 1 stellar mass-grid, we model a typical CoRoT- 
field with a mixture of 56 %, 37 %, 7% of 2, 1, and 0.8 M Q 
stars, resp.. This roughly represents a grouping o f stellar spec- 
tral types A + F, F + G and K, resp.. As noted bv lBroegl {2006) 
Mark 1 has a relatively sparse sampling of the stellar mass- 
dependence but is already based on ss 10 6 planetary models und 
thus by far the largest survey of planet-formation. We do not in- 
clude the Mark 1 results for 0.4 Mq (M-stars) here. To account 
for the uncertainties in planetesimal-formation and -accretion 
theory we use the results for all core-accretion rates given in 
Mark 1, i.e. 10~ 6 , 10~ 4 , and lO^M^a- 1 . The resulting the- 
oretical planetary mass spectra for a CoRoT-stellar population 
mix, the CoRoT-Mix for short, are shown in Fig.Q] We note that 
this a-priori prediction of the planetary initial mass function nat- 
urally predicts known empirical 'planetary' masses. Given the 
fact that the only particular macroscopic input masses are the 
stellar masses, that is non-trivial. Fo r a general discussion of the 
properties we refer to Broeg (2006) and focus on the particular- 
ities of the CoRoT-Mix in the following. 

The multi-peak structure in the mass fun ctions follow s the 
general double peaked structure described by Broeg (2006) but 
is complicated by mixing the mass spectra of different host stars. 
This results is the structures displayed in Fig. Q] Averaging the 
mass spectra using a finer sampling in host star mass would most 
likely revert the mass spectra to a double peaked structure. At all 
periods the most frequent planets appear with the first peak near 
Jupiter's mass and the second one close to Neptune's. Both peaks 
move to larger masses with increasing period or orbital distance. 
The shift is about half an order of magnitude. The peaks in the 
mass spectra become narrower for larger distance indicating that 
the typical planetary masses (Neptune and Jupiter) become more 
well defined farther from the host stars. A key feature at all pe- 
riods is a rich population of 'Hot Neptunes' that always signifi- 
cantly outnumbers the Jupiters. 
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a Estimated from lBroed j2006l) . Figs. 8.2, B.3 and B. 1 1. 



Table 1. Characteristic quantities of the CoRoT-Mark 1 vl.l 
theore t ical pl anetary population for a solar mass host star, from 
iBroegl (120061) for various orbital periods, r o ,b. Semi-mayor axis, 
a and the radiative equilibrium nebula-temperature, T ne b result 
from the host star mass and luminosity, resp.. M\ and M2 give 
characteristic values for peaks in the planetary mass spectra. Pi 
and P2 are the corresponding typical nebula pressures. 



For reference and an accompanying study of the impact of 
particle mass-loss processes characteristic numbers are given for 
the pure solar mass results of CoRoT-Mark 1 in Tab. 1 . 

4. Stability of the planetary population 

The CoRoT Mark 1 planet population is based on an equilibrium 
concept, much like the stellar main sequence where stars are in 
complete i.e. hydrostatic and thermal equilibrium. The advan- 
tage is that whatever the formation history of the planet and the 
nebula was, it does not influence the planetary properties. Thus 
whether planets migrated to their final position or had undergone 
giant impacts or mass-loss events during the formation era does 
not change our prediction as long as the respective conditions 
of the planetary equilibria are realised — i.e. nebulae are suffi- 
ciently diverse — as the planets approach their final masses. The 
disadvantages are (1) that because the method is based on time- 
independent models, we do not automatically have information 
about the stability of the equilibria used in the construction and 
(2) because it is based on planet formation theory we do not ac- 
count for modifications of the planetary properties during plane- 
tary evolution in to the present, in particular and most important 
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mass loss, that might distort the mass spectra. The second issue 
is discussed in an accompanying letter, the q uestion of stability 
in the following sections. For HD 149 026 b. iBroeg & Wuchterll 
(2006) showed that planet formation is a completely hydrostatic 
process in that case. They gave a complete formation-history in- 
cluding a radiation fluid-dynamical calculation with convection 
and detailed, particle in box planetesimal accretion to ages of 
about 100 Ma. Because that is presently beyond computational 
feasibility for the entire Mark 1 population of ~ 10 6 planets 
we have addressed the stability question by a series of simpli- 
fied model systems: (1) linear analysis of isothermal and poly- 
tropic protoplanets, (2) non-linear fluid dynamics for isothermal 
protoplanets, (3) quasi-hydrostatic stellar evolution type calcu- 
lations to investigate non-linear secular or thermal stability and 
(4) radiation-fluid dynamical calculations with convection. 



5. Connections between mass spectra and linear 
dynamical stability analysis 

A compre hensive line a r analy sis of isothermal protoplanets was 
given by Sch onke.TI (|2005) who found two types of linear- 
instabilities that operate in protoplanets for sufficiently high- 
nebula-densities and are damped by large cores. One charac- 
teristic peak in the mass spectra obviously results from struc- 
tures in the outer envelope, which are in a state describable as 
the onset of self-gravity. These models were found for a cer- 
tain region in the parameter set of core mass and gas pressure at 
the core surfa ce (cf. 'region IV in IPecnik & Wuchterll d2005l) or 
Broeg (2006)) and have all the same characteristic mass, lead- 
ing to the peak in the spectra. The linear stability analysis (cf. 
Scho nke.TI (|2005)) shows a clear connection between this class 
of models and a dynamical instability. Almost all models are 
linearly unstable, except for a small strip in the low pressure 
regime. What kind of consequences for the spectra follow from 
this result can not be answered completely from the standpoint 
of the linear analysis because no information about the dynam- 
ical evolution can be gathered from it. It is likely that a part 
of the unstable models will end up in the stable low pressure 
strip, a behaviour which the non-linear stab ility analysis shows 
(cf. lPecnik! d2005l) . iPecnik & Wuchteril d2006l) ). At least we can 
say that a 'stability-corrected' mass spectrum would reduce the 
'strength' of the peak significantly. The appearance of the de- 
scribed unstable models mainly depends on the available enve- 
lope volume (i.e. the Hill radius) and therefore on the orbital 
distance or period. As the Hill radius becomes small (i.e. orbital 
distance or period are small), the region where unstable enve- 
lope structures are built up shifts to smaller core masses and the 
spectral peak broadens. Because of the quantitative differences 
between isothermal and detailed non-isothermal protoplanets it 
is nor possible to presently decide which peak in the Mark 1 
spectra will suffer the flattening, likely the more massive one of 
the two main peaks but possibly both. 



6. Non-linear stability of short-period isothermal 
proto-planets 

Isothermal solution manifolds Pecnik & Wuchterll d2005l) for 
short-period (P=1..4d) proto-planets (PPs) have different non- 
linear stability properties than tho se of long-period PPs (see 
IPecnik. B. and Wuchterl. G] d2007l) for dynamical properties 
summary, and Pecnik (2005) for a — 5 AU orbital region de- 
tails). 



Unlike long-period PPs, the critical core mass (M cl ; t ) for 
short-period isothermal PPs gets larger for smaller orbital pe- 
riods. Isothermal Hot-Jupiter's M a i t also strongly depends on 
the temperature and the density of the surrounding nebula. 
Additionally, massive subcritical short-peri od PPs are more lin- 
early stable than long-period counterparts lSchonke. J.I d2005l) . 

Fluid-dynamical calculations at orbital periods of 1 and 4d 
for a solar mass primary, show that the non-linear stability 
properties of isothermal subcritical PPs with temperatures of 
T\ = 1260 K and 72=1600 K are almost identical, which hints that 
most of these equilibria might be thermally stable. Short-period 
subcritical PPs with gas density at the core surface (g cs ) of up to 
around 10 kg irT 3 are stable and just oscillate around the equi- 
librium if perturbed. Subcritical PPs with higher g cs will make a 
transition towards a corresponding solution with lower g cs , and 
could either loose or gain envelope mass in the process. Solution 
with still higher p cs will eject more than 95% of its envelope. 
Finally, a compact PP (a PP with highest g cs ) will remain stable, 
but will include an outgoing wind, supersonic at the outer bound- 
ary. Depending on the level of envelope compactness, mass loss 
due to this wind could be anything from completely negligible to 
a significant percent of the envelope mass within a single sound- 
crossing time. Additional effect of this wind is "puffing-up" of 
the outer PP's stratifications. All of these properties suggest that 
massive stable subcritical PPs are formed more easily at small 
orbital distances (than at larger orbital distances), and that quasi- 
hydrostatic evolutionary paths are available up to the M cl -i t . 

7. Quasi-hydrostatic calculations and 
secular/thermal stability 

The previous sections were dedicated to the isothermal case that 
is computationally very efficient but cannot provide quantitative 
agreement with present day giant planet masses. To analyze the 
stability of the Mark 1 population directly, i.e. with the time- 
dependent versions corresponding to the the static Mark 1 equa- 
tions, a full stellar evolution type calculation is necessary. Such 
calculations, with constitutive relations identical to the Mark 1, 
have been performed for selected parts of the planetary mani- 
folds for orbital periods of 1 and 4 days for the solar mass case. 
The result is that in all cases calculated, i.e. for core masses to 
10 23 kg and total masses up to about 50 earth masses the pro- 
toplanets contract and loose less than 30% of their total mass 
within the 100 Ma calculation period. 

8. Non-linear convective radiation-fluid-dynamic 
calculations 

Finally we performed full nonlinear-radiation fluid calcula- 
tions for selected proto planetary models in continuation of 
Broeg & Wuchterl (2006). These calculations are the largest 
computational effort and have been done in parts of the CoRoT- 
Mark 1 planet population that are thought t o be p a rticula rly un- 
stable from both the population analysis, Broeg (120061) which 
finds the planets at largest orbital distance to more likely become 
dynamical a nd th e isothermal li near and non-linear studies, 
ISchonke. J.I d2005l) : IPecnik! (120051) . that find the high envelope 
mass, high nebula pressure, low core mass protoplanets to be the 
most likely ones to suffer an instability. Because we are pr i marily 
interested in the stability here, unlike Broeg & Wuchterll ([2006) 
we use a constant nebula state for the outer boundary condition 
— in practice that means an unlimited mass reservoir instead of 
a feeding zone with given mass. 



4 



G. Wuchterl et al.: Properties of the short period CoRoT-planet population I: 



We used the equations of IWucht erl & Tschar nuterl d2003l) 
with planetary boundary conditions, IWuchterJ d 199 lb . The dy- 
namic evolution of 6 planets of the Mark 1 population for a 16d 
period at 1 Mq and core masses of 0.47 and 1.11 M eart h was cal- 
culated with static models as initial condition. These core masses 
are relatively low c ompared to th e typical critical masses for 
close-in planets, see iBroegl (120061) to select potentially unstable 
planets. It turned out that the envelope mass average of the first 
adiabatic exponent, F\, an indicator for dynamical stability in a 
stellar structure type context, has a characteristic depletion be- 
low the critical value of 4/3. T hat occurs as th e 'island' of large 
planetary envelope masses (cf. Broeg (2006)) is crossed along 
constant core mass. 

In all six cases the evolution remains quasi-hydrostatic and 
the masses within 10% of the initial, static value for 2000a. 
Dynamically this means that all the respective protoplanets are 
stable. That in spite of selecting two models by looking for the 
lowest Fi at given core-mass. There is also no indication of rele- 
vant mass loss due to a secular instability. The two most massive 
planets with 105 and 126 earth masses, initially show a few per- 
cent reduction in mass as the initially static planet starts to con- 
tract. After a thermal adjustment of a 10 4 a the mass levels off at 
a few percent less than the initial values and starts growing. 



9. Close-in Convective Protoplanets 

The fact that 'realistic close-in exoplanets turn out to be more 
stable than estimated by isothermal models can be explained 
the fact that the env elopes are largel y convective . Largel y con- 
vective protoplanets. lWuchterj dl993l) : llkoma et alj d200ll) damp 
dynam ical instabilities in their lar ge convection zones, IWuchterll 
(1995) and are favour accretion. IWuchterl et alj d2000t) . 

Thus as far as the sample calculations are representative of 
the ensemble, the planets of the CoRoT Mark 1 population are 
stable. The onset of apparent adiabatic dynamical instability in 
2/6 models is damped by convection. 



10. Discussion 

The CoRoT-Mark 1 theoretical planet populations rely on the 
implicit assumption that there is no preferred scale in the masses 
of planetary cores and no preferred range of nebula pressures. 
Together this requires to assume a diversity of gravitationally 
stable nebulae. That does not imply arbitrary ranges for core- 
masses and nebula pressures but sufficiently wide ranges so that 
the planetary core-envelope equilibria can 'draw' any required 
value from a nebula at a given orbital radius. If some cores that 
are required in the population do not form, the respective plan- 
ets would not be able to form and thus not appear in the ob- 
served mass spectra. Presently we do neither have a sufficiently 
comprehensive overview about the diversity of nebula proper- 
ties from observations nor from theory. Thus the CoRoT-Mark 1 
predictions rely on an assumption of nebula diversity. Presently 
the best argument that there is indeed a huge diversity of nebular 
conditions comes from the observation of extreme planets as GQ 
Lupi b and HD 149 026 b. 

The second key assumption to derive the mass spectra is that 
the planets are static. While planet formation has been found to 
be a hydrostatic process under many conditions and for most of 
the time, there are important exceptions and there is no guarantee 
that the results from studies of planet formation for a rather re- 
stricted set of conditions can be generalised to the diverse nebula 
properties required by the Mark 1 population. 



A general study of stability for the entire ensemble of ~ 10 6 
'realistic' protoplanets is presently beyond our computational 
reach. We therefore followed a two-fold strategy. We studied 
simplified problems (linear isothermal stability) in a compre- 
hensive way and sampled a hierarchy of more detailed models - 
nonlinear dynamical stability of isothermal PP, nonlinear secular 
stability of realistic PP and did a limited set of non-linear radia- 
tion fluid-dynamical calculations at positions that were indicated 
as particularly unstable in the more comprehensive but simpler 
models. While the simpler models provide a detailed and synop- 
tic view, they indicate instabilities in some parts of the planet- 
manifolds. Approximately 1/3 to 1/6, say, of the more massive 
envelopes with small core are unstable and likely correspond- 
ing to realistic Mark 1 planets. The few detailed models show 
that improved physics, in particular time-dependent convection 
acts stabilising at the investigated places an indicates stability 
for all cases calculated. In particular a non-linear development 
of a case with linear adiabatic dynamic instability indicated by 
globally low first adiabatic exponents (down to 1 .25) is damped 
out by convection and the respective planets remain static and 
thus stable, valid elements of the Mark 1 planet population. 

11. Conclusions: CoRoT's planets 

Based on the assumption that diverse gravitationally stable pro- 
toplanetary nebulae are produced by the star formation process 
we predict the short period, r or b < 16d planet population of a 
typical CoRoT field: (1) they closely follow the CoRoT-Mark 1 
planet population applied to the stellar population mix of the 
field, i.e. show planetary initial mass functions (IMFs) as given 
in Fig. [TJ (2) these IMFs have to be corrected for particle loss 
process to arrive at the epoch of observation mass function; this 
corrections are discussed in a parallel letter; (3) typical field 
planetary IMFs depend on orbital period an d have t he la rgest 
mass peak around a Jupiter mass, we refer to iBroed d2006l) for 
more predictions of properties of single-host-star-mass planetary 
IMFs. (4) the broad peaks in the stellar population mix IMFs 
point to a large population with masses of 1 - 3 10 26 kg, (50 - 
160 earth-masses). (5) there is an abundant population of Hot 
Neptunes and an important population filling the solar-system 
mass gap between Neptune and Saturn. 
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